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We consider in detail the -dependence of the DIS structure functions, with Q being the 
virtual photon momentum. Very often this dependence is claimed to be originated by the Q^- 
dependence of the QCD coupling. This leads to the small- a; asymptotics of the structure functions 
with -dependent intercepts. We demonstrate that the DGLAP parametrization ots — as{Q^) 
is an approximation valid in the region of large x (where 2pq can be approximated by Q^) only, 
providing the factorization scale is also large. Outside this region, the DGLAP parametrization 
fails, so Os should be replaced by an effective coupling which is independent of at small x. As a 
consequence, intercepts of the structure functions are independent of Q^. Nevertheless, the small- a; 
asymptotics of the structure functions explicitly depend on Q^, even when the coupling does not 
depend on it. We also consider the structure functions at small and give a comment on power-Q^ 
corrections to the structure functions at large and small Q^. 
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I. INTRODUCTION 

Each of the DIS structure functions has two independent arguments: the invariant energy w — 2pq (with p and 
q being the hadron target and virtual photon momenta respectively) and virtuality of the incoming photon = 
—q^ > 0. Early experiments of 60s at SLAC proved that the structure functions depended on only one argument 
X — Q'^/w. This phenomenon was called scaling. Later, new experiments proved that the scaling was violated i.e. in 
addition to a;-dependence the DIS structure functions also depended on straightforwardly. Since that it become 
conventional to regard variables x and as independent arguments of the structure functions. However, one has to 
realize that those arguments (in contrast to the set of really independent arguments w, Q^) cannot be regarded as 
independent variables at all times. Indeed, they are independent when is fixed while w is scanned but they are 
not independent at fixed w and scanned Q^. The -dependence of the structure functions is the objective of the 
DGLAP evolution equations [l| implementing the total resummation of In'^ Q^. In DGLAP and beyond this approach, 
the -dependence of the structure functions is often thought to be induced by the -dependence of as'- 

a,^a,{Q^). (1) 

As this parametrization takes place in the DGLAP equationspj], throughout the present paper we will address it 
as the DGLAP-parametrization. Parametrization ([IJ, originally deduced for large x, where 2pq can be approximated 
by Q^, is often exploited to describe DIS in the small- a; region and even in the small- a; asymptotic of the structure 
functions, i.e. at a; — t- (see e.g. Refs. As a result, the small- a; asymptotics look like 

f-^-^^'''\ (2) 

with the Q^-dependence of the asymptotics originated by the "intercept" A(Q^). 

In the present paper we argue against both validity of Eq. ([2]) and treating Eq. ([1]) as an exact formula which 
could be used at any x and Q^. In what follows we will show that the parametrization ([1]) should be used with a 
certain care: it can be used within the proper DGLAP kinematic region, i.e. at large x and Q^, however it is not 
enough: the factorization scale should also be large. We show that if any of these requirements is violated, asiQ"^) 
should be replaced by the effective coupling which does not depend on at small x as well as at small Q^. The 
effective coupling, being independent of at small x, cannot bring any Q^-dependence to the structure functions. 
In this case the -dependence of the structure functions is originated by alternative sources. We also argue against 
regarding the approximative parametrization ([l} as the basis for the shift as{Q^) — >■ cXs{Q^ + M^) ^-s a continuation 
into the small-Q^ region when the shift is used in order to keep as in the perturbative domain. Instead we suggest 
an alternative approach to describe the structure functions in the small-Q^ region. 

Our paper is organized as follows: In Sect. II we consider the DIS structure functions in general, focusing on 
their conventional dependence on Q^. In Sect. Ill we consider the structure functions in the framework of DGLAP 
and discuss parameterizations of the QCD coupling in different kinematic regions. In Sect. IV we derive a general 
expressions for the structure functions at small x and derive a general form of their small-x asymptotics, demonstrating 
that there is no room for Eq. ^ a.t x 0. In Sect. V we discuss the role of singular factors conventionally used in 
DGLAP-fits for initial parton distributions. We consider the structure functions at small in Sect. VI. We make a 
brief remark on higher twists in Sect. VII. Finally, Sect. VIII is for our concluding remarks. 

II. GENERAL STRUCTURE OF THE DIS STRUCTURE FUNCTIONS 

All DIS structure functions are given by a convolution of perturbative components J^^'^'^*-' and parton distributions 
^r- For instance, in collinear factorization and when DGLAP is used, the non-singlet structure function fNs{x,Q'^) 
is represented as the convolution 

fNs{x,Q^)^Cy{x,Xo)®Aq{xo,Q^) (3) 

where Cq is the quark coefficient function while Aq{Q^) stands for the quark evolved distribution defined at scale Q^. 
Cg(x,xo) controls evolution of this distribution from xq to x. The -evolution of the quark distribution from the 
initial scale /i^ to scale is controlled by the DGLAP evolution equations [l]. It can be written as 



(4) 
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where the initial quark distribution Sq(xo,is^) is conventionahy chosen at ccq 1 and fi ^ 1 GeV. The scale fi^ is 
called the factorization scale. Eqs. ([3J U]) can be written as one convolution: 

where /jfg'^*^ denotes the whole perturbative evolution of 5q both with respect to x and Q^. Eq. ([S|) can easily be gen- 
eralized to the case of the singlet structure functions by adding gluon contributions. So, both in basic factorization^ 
as well as in collinear and -factorizations, such a convolution can be written as 

/(x, = ^ /t^'-*) (q^, Q\ K^) ® {pn, K^) , (6) 

r 

where $ denotes parton distributions in any of the factorizations, the subscript r marks the intermediate partons 
(quarks and qluons) with momentum k and the summation over r takes place when necessary. The notation ® 
means the one-dimensional integration over the longitudinal for collinear factorization and two-dimentional (both the 
longitudinal and transverse momenta) in the case of -factorization; = ^j? in the case of collinear factorization. 
In Eq. ([B]) we have dropped a possible dependence on unessential variables like the target spin and mass. 

A. General structure of the perturbative component 

The simplest form of structure functions corresponds to collinear factorization where it can be written in terms the 
Mellin integral: 

, r°° du: _ 
f = X 



— too 



flf-^'\u:)5q{u:) + /(f^*' (w)^3(^), J (7) 

with 5g and 5g(Lo) standing for the quark and gluon distributions in the uj -space respectively. In Eq. ([7]) they are 
convoluted with the perturbari^ 
the following generic structure: 



convoluted with the perturbarive contributions fr^'^^*'\ These perturbative contributions in any factorization have 



/200 T _ 
— :E-"a(c^, A^2)e"^(-.«-Q^«') (8) 
-loo 27ri 

where the term Cr is a generic notation for the coefficient functions and Sl^ is expressed through the matrix of 
anomalous dimensions. For instance in LO DGLAP, where, in collinear factorization, — /i^, the non-singlet 
perturbative component fj^g^^\= /q^'^'^*') is expressed in terms of the non-singlet coefficient function Cms and ^Ins- 



C 



In (QVA^) 
In (m7A2) 



(9) 



with jq'l\uj) being the well-known quark-quark anomalous dimension and b = (33 — 2n/)/(127r). In NLO DGLAP, 
both Cns and fl^s acquire additional terms proportional to a^. In particular, the structure of the NLO coefficient 
functions is 

c'^Ns'^'''^ = l + a,iQ^)CNsi^)- (10) 

We have used this simple example in order to demonstrate explicitly that the -dependence of any structure 
function / is achieved through the -dependence of its perturbative components fr^'^^^\ The -dependence of 
jip&rt) j,pj-|^gg fj.Q]-|2 ^j^g upper limit of integrations over fc^ in expressions for both $7^, as explicitly shown in Eq. ([9]), 
and in expressions for coefficient functions. The latter is related to the parametrization as = as{Q'^). This feature is 
common for all structure functions, although specific expressions for C^r' and firr' are different for different structure 
functions. 
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III. q2 -dependence of the qcd coupling 



In this Sect, we discuss parameterizations of in DGLAP equations in collinear factorization at different values 
of X and fj? Although for the sake of simplicity we consider the LO DGLAP equations, our conclusions are valid for 
NLO, NNLO and so on. In the LO DGLAP evolution equation the parametrization ([1]) appears when the integral 
DGLAP equations for /(p'^''*) 

/(p-*)(a:,Q2)^ ^ J' ^aikl)f^r'\x/P,Qykl)PMP) (H) 
are reduced to the differential equations 



Eqs. (|llll2p demonstrates explicitly that Eq. ([T]) holds when 

(i) the upper limit of integration over fc^ in Eqs. (jlip is which corresponds to x 

(ii) the parametrization 



as = as(fci) (13) 

is used in the involved Feynman graphs. Throughout the present paper we will refer to the parametrization of 
Eq. (fT3|) as the standard parametrization. According to the results of Ref. [a| , one can use the standard parametrization 
([T^ in expressions for j'(P'=''*) in collinear factorization only when the following two conditions are fulfilled: First, values 
of X should be large: 

X ~ 1 (14) 

to ensure the use of instead of 2pq. Second, the factorization point should obey the strong inequality 

fi^ > e'^A^ w 23A2. (15) 

When one of both of the requirements in Eqs. [T5|) are violated, the coupling as(A;^) in the integral DGLAP 
evolutions equations for f(p<^^*)^ should be replaced by the effective coupling ae//, as shown in Ref. [1]. This coupling 
can be represented by different approximative expressions, depending on the kinematics. In Subsections A and B we 
consider the cases of violating either or (ITSI) and consider violation of both of them in Subsection C. 



A. Case A 

When Eq. (jTSl) is violated but x is so close to 1 that essential /? in Eq. pT|) are ~ 1, the coupling as{k'^) in Eq. pT|) 
should be replaced by 

where we have denoted 

Zo = hV/A'), r(fci) -ln(fci/A2). (17) 
It leads to replacement of the integration of ag (fc^) in Eq. (|9]) by the integration of «£/ f{{k±) over the same interval: 

Let us notice that when fj, obeys Eq. (|15p . the first and second terms in Eq. PB|) cancel each other while the third 
term can be approximated by as{k'^). This reduces the integral in Eq. (IT51) to the DGLAP expression in Eq. ([5]). 
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B. Case B 

Let us consider the opposite situation when /i satisfies Eq. (flS)) but x is small. In this case as(fc^) in Eq. (fTTj) should 
be replaced by Q;s(fc^//3) and the upper limit of integration over fc^ should be changed. The point is that integration 
over fc^ always runs from some starting point fj? (which is fixed from physical considerations and sometimes coincides 
with the factorization scale) to the total invariant energy 

s^{p + qf ^w{l-x), (19) 

with w — 2pq. In the DGLAP framework, where x is not far from 1, s ?s Q'^, so in the DGLAP equation (|lll) the 
upper limit is Q^. On the other hand, s « w at small x. This leads to replacement of Eq. (|11|) by the following 
equation: 

/(p-*)(^,q2) ^ r r ^o.s{ki/p)fir'\xip,Q^/ki)PMP). (20) 

with 

Po^x + ki/w. (21) 

Obviously, /?o ~ a; at large x. In addition, one can neglect the j3 -dependence of as{k\/ P) at large x, arriving back 
to the integral DGLAP equation (fTTj). 

C. Case C 

Finally, when ^ does not satisfy Eq. p5|) and, in addition, x <C 1, the coupling as{k']_l (i) in Eq. should be 
replaced by aef f{k\/ P): 

aeff{kl/P) = \jfT^^ - ^ arctan + ^ arctan (^) , (22) 

with I — In (fc^/(/3A^)). It converts Eq. (|20| into a new equation: 

/(^^'■*)(x,Q^) ^ T'^f fas{kl/p)flr'\x/P,Q'/kl)PMP). (23) 

Eqs. ()20l23p explicitly demonstrate that there is no factorization between the integrations over (3 and k'j^ at small 
X. In other words, the factorization between the longitudinal and transverse spaces taking place in DGLAP vanish 
in the small- a; kinematics. Eqs. ()20|23p also show that structure functions at small x depend on through the 
integration limit /3o while the effective couplings at x ^ 1 differ a lot from the standard parametrization given by 
Eq. ([T|). Eq. corresponds to the case when both Eq. (|14l) and are violated while the expressions for the 
couplings and structure functions in Subsects. A and B correspond to violation of either Eq. (IT4t or (fT5|) and can 
easily be obtained from Eqs. (|22I23|) . Nevertheless, we accentuate that Eq. is also approximation obtained in 
Ref. in order to factorize the coupling. This expression should be used only in the context of evolution equations. 
A more general treatment of the QCD coupling, where the factorization was not required, was done in Ref. [9]. We 
discuss it briefiy in the next Sect. To conclude, we would like to notice that Eq. ((23)) was obtained with constructing 
Dispersion Relations for forward scattering amplitudes. This approach is similar to the approach of Ref. [l3| but 
differs from the alternative approaches (see e.g. Refs. [Ill where the Dispersion Relations were constructed for 
as by itself. 

Obviously, logarithms of x are large in the small-x region, so their resummation to all orders in as is important. 
This is beyond of the reach of DGLAP, where only logarithms of are resummed. A generalization of the DGLAP 
equations to the small- a: region was done in Refs. JJJ] (see also overview JAl) for the spin structure function 171 and 
non-singlet component of Fi : there were resummed the leading, double logarithms of both x and and at the same 
time the QCD coupling in each rung of the involved Feynman graphs was running. Applying the saddle-point method 
to the structure functions obtained in [Tsj exhibited their Regge behavior at a; — > 0, with the intercepts being just 
numbers without any parameters (e.g as, etc). The intercept of the non-singlet Fi proved to be 0.38, the intercept of 
the non-singlet gi was 0.42 and the intercept of the singlet gi was 0.86. 
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IV. GENERAL STRUCTURE OF THE EVOLUTION EQUATIONS FOR DIS STRUCTURE 

FUNCTIONS 

In this Sect, we, consider general expressions for the structure functions, skipping unessential details. Throughout 
the present Sect, we imply that coUinear factorization is used, though a generalization to kx -factorization is easy to 
do. As equations for the structure functions involve convolutions (see e.g. Eq. (O), it is convenient to represent the 
perturbative contribution fl^'^'^^'^ in terms of the Mellin transform: 

^x-Fr{co,y), (24) 

were we have denoted y = ln((5^//i^) and introduced the Mellin amplitude F^. Focusing on the t -channel intermediate 
state with two partons and neglecting other contributions (see for detail Ref. [14|) allows us to compose the following 
system of equations for F^: 



[uj + d/dy]Fr{oj,y) = Fr'{uj,y)Hr'r{oj) 



(25) 



Solving these equations, we arrive at 



^x--Fr{oj, y)Sr{u;) = ^ / |^^x-a(^)e^"'-(-)fr(c.) 



(26) 



where r = q,g, so that 6r = Sq,dg, and ^Ir is made of new matrix of anomalous dimensions i/,./,.. Both Cr and Hr>r 
account for both the total resummation of logarithms of x and running as effects. As as depends on the longitudinal 
Sudakov variables, it participate in the Mellin transform (see Ref. Q for detail). As a result, as with the time- like 
argument is replaced by 



7/^ 



dpe 



(27) 



where 77 = \a{ii^ / K'q(j j^) and ji is the IR cut-off. When the argument of the coupling is space-like, the tt^ -terms are 
absent. 

Comparison of Eqs. (|26p and ^ shows that the coefficient functions Cr do not depend on at all. The 
-dependence of fr{x,Q^) is located in the exponent of (|26p and it is controlled by the anomalous dimensions. The 
small-x asymptotics of fr can be obtained with applying the saddle-point method to Eq. (j26p . The stationary phase 
is obtained as a solution to the following equation : 



— K + ina(c^) + yf^r(c^)] = e + 7^-7^ + y-T^ = 0, 

dw Cr du! dcu 



(28) 



where ^ = ln(l/a;). This equation should be solved at ^ 00 and fixed y and the solution, wq, is called the stationary 
point. The rightmost root of Eq. stipulates the small- a; asymptotic behavior of fr, with 



A — niax[a;o] 



(29) 



being the intercept of the structure function. Eq. (1^51) implies that the large factor ^ must be equated by another 
large factor. As y is fixed, such a factor can be any of the following options : 



(i) 

(ii) 
(iii) 



Cr{uj) -0, 

dCr{uj)/duj — > —00, 
dQri'-^) / dio —00, 



(30) 
(31) 
(32) 
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when Lu ^ uiQ. Accounting for logarithmic contributions ~ In" a; to Cr{uj) and ilr means that they acquire contribu- 
tions ~ c„/ti;^+" and ~ c^/w" respectively and therefore any of them and dCr/dui, which is ~ c„/w^"'"", are singular 
a.t u! —i' ujQ — 0. Obviously, fulfilment of Eq. (pO)) can be achieved only when a part of Ck is negative, which contradicts 
to calculations made in DGLAP and to expressions for the non-singlet structure functions obtained in Ref. [l^. Then, 
let us note that as contributions to dCr/duj are more singular than contributions to Cr or fir, the case (ii) is the 
most important compared to (i) and (iii). Indeed, ut should be small in order to prevent oscillations of the factor 
exp(^CL') in Eq. (I26|) at large ^. The series of the pole contributions ^ 1/lu'' can be summed up. For example, the 
coefficient function Cns a-nd anomalous dimension Hns for the non-singlet structure function F{^^ (for this case 
flr{uj) is replaced by HNsi'-^)) proved to be 



Cns 



oj - Hns ' 



Hns = (1/2) U - - ^a'5(^) 



(33) 



When tts is fixed, Bns — ^asCp/n, otherwise it is given by a much more complicated expression (see Ref. jl3| for 
detail). It is easy to see that all coefficients in the expansion of Cns into series in l/w" are positive. This also excludes 
Eq. ([30]) and demonstrates that Eq. (|3T|) corresponds to the most important case. So, we arrive at the asymptotics of 
/ in the form of following contributions of the Regge type: at x — > 



fix,Q^)^x-^{Qyf,^f^^Ur. (34) 

Eqs. (|26I34I) demonstrate explicitly that the dependence of the structure functions on has nothing to do with 
the -dependence of the coupling. 



V. REMARK ON THE DGLAP-FITS FOR STRUCTURE FUNCTIONS 



The asymptotics ([34]) is of the Regge type. The Regge behavior is generated by total tesummations of leading 
logarithms of x and it is unrelated to behavior of the initial parton densities Sq{x) and Sg{x). In particular, both 
Sq{x) and Sg{x) are not supposed to reveal a power behavior ~ x^" at small x. On the contrary, the standard DGLAP 
fits 

Sq,Sg = Nx-''{l~x)''{l + cxy, (35) 

with a,b,c,d,N > 0, conventionally includes such factors. The singular factors x^" are incorporated in the fits ad 
hoc in order to provide / with a steep growth at small x and thereby meet experimental data. Indeed, in the lo -space 
the factor corresponds to the pole contribution 

= l/{uj — a). (36) 
Being installed in Eq. ([7]), the pole contribution ([55]) straightforwardly leads to the Regge asymptotics 



f^x-"" 



ln(gVA2) 



7(a)/6 



ln(/i2/A2) 

which differs a lot from both Eq. (|34| and the well-known LO DGLAP -asymptotics 



(37) 



'■■»/-)'" [Iffy] ^ P«) 

Let us note that Eqs. (I34I38P involve non-singular fits. Confronting to ([57)1 proves that the only role of the 
singular factors x~°' is, actually, to mimic the total resummation of In" x. When the resummation is accounted for, 
these factors can be dropped and therefore the number of parameters in the fits can be reduced. The qualitative 
difference between the asymptotics (p4)) and ([37| includes two items: 

(i) The intercept A in Eq. (|34l) is calculated while a in p7l) is fixed from experiment. 

(ii) The factor x~'^ in Eq. (IMl) is present in the asymptotic expressions, i.e. at a; — 0, only. It never appears in 
expressions for the structure functions like Eq. ([T]) at finite x. On the contrary as soon as the fit (1351) is used, the 
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factor x~°' is present in Eq. (jT]) at any cc, including large x ^ 1. It means that, in a sense, the x -dependence of 
/ is always controlled by its asymptotics, when singular fits are used. This contradicts to the observation made in 
Ref. [13]: the small- a; asymptotics represent the structure functions reliable at very small x [x < 10^^) only while at 
X > 10~* the structure function is much greater than its asymptotics. The latter leads to the quantitative difference 
between a and A: In order to force Eq. (|37p to represent / at presently available x, i.e. at a; > 10~^, one has to 
increase a, which leads to the following relation: 

a > A. (39) 

By this reason, we will name the fictitious intercepts a as pseudo-intercepts. As their values are determined from 
experiment at available x, they always exceed genuine intercepts. So, the situation looks like that: On one hand, 
the use of the pseudo- intercepts a allows one to approximate Eq. ([7]) by its asymptotics which is given by simple 
Regge-like expressions ([37]). On the other hand, the power x -dependence in Eq. (f37|) has nothing to do with impact of 
genuine Reggeons introduced in Theory of the Regge poles. The wide-spread tactics to use pseudo- Reggeons, instead 
of resummation of Inx, in order to solve immediate practical tasks leads to serious theoretical problems especially 
important for the singlet structure function Fi and the singlet parton distributions: there are no available expressions 
for those objects besides their Regge asymptotics x^^'' , with Ap being the Pomeron intercept. In order to fit such 
asymptotics to explanation of experimental data at available energies, they use one or several Pomerons with large 
pseudo- intercepts Ap violating the Froissart bound. Finally, let us note that necessity to include singular factors in 
the fits is a clear indication that essential logarithms of x are not accounted for. 

VI. STRUCTURE FUNCTIONS AT SMALL 

Description of the structure functions at small is important because this kinematics has been investigated ex- 
perimentally. For instance, the spin-dependent structure functions at small are investigated by the COMPASS 
Collaboration (see e.g. Ref. [IBl)- On the other hand, this region is absolutely beyond the reach of DGLAP. De- 
spite this, sometimes in the literature (see e.g. Refs. 0, recent paper Ref. [![ and refs therein) the DGLAP 
-parametrization — as{Q'^) is treated as an exact expression where can acquire an arbitrarily small values and 
therefore the shift 

a,(g2)^a,(g2+^2) 

is needed at small and especially at — )■ in order to keep as within the perturbative domain. However, in 
the present paper we have shown that the DGLAP-parametrization = as((5^) is an approximation valid for large 
only, where the shift (PH)) is totally unnecessary and even cannot be seen at ^j? <ti ■ Actually, the DGLAP 
-parametrization fails at small or at small x , so as{Q^) should be replaced by the effective coupling Ueff- The 
coupling ttef f does not depend on at small or small x at all, which makes the shift (|40p absolutely unnecessary. 
Let us stress that interpretation of the approximation Us — as{Q^) and the shift (^0]) as exact expressions has led to 
various misconceptions abundant in the literature. On the other hand, the small-Q^ kinematics has been investigated 
experimentally by the COMPASS Collaboration (see e.g. Ref. [11]), so it is important to describe the structure 
functions in this region. In Ref. [l^ we proved that the expressions in Ref. [l^ for the structure function gi at small 
X and large can be extended to small by shifting 

kl ^ kl + ^? (41) 

in propagators of the soft quarks and gluons. The shift (|41l) . where is a cut-off, is necessary to regulate infrared 
singularities in involved Feynman graphs. This shift eventually leads to the shift 

Q2 ^ Q2 ^ (42) 

in logarithmic contributions to the structure functions. Our estimate for ^ obtained with using the Principle of 
Minimal Sensitivity [3| was /i « IOAqcd (for more details see our recent paper [19]). The shift (|42]) can be neglected 
at » fj? but becomes essential at small Q^, making possible description of the kinematic region where are 
small. For example, DIS at such a kinematics has been studied experimentally by the COMPASS collaboration. 
Therefore, it is the shift (HIT) that brings theoretical grounds for description of the small -region. Let us stress that 
this shift should not been applied to as and it has nothing to do with the baseless modification (|40)) of the coupling. 
On the other hand, the shift (|^T|l together with the shift Usik^) — > ^^(fc^ -I- IJ?) can be used when k are momenta 
of the virtual t-channel gluons connecting the perturbative part of the structure function to the parton distributions 
but these shifts have nothing to do with the shift of Eq. (PH)) (see Ref. [l^ for details). 
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VII. REMARK ON HIGHER TWISTS 

It is interesting to notice that the shift (|40)) can be used to clarify the problem of higher twists and, in addition, it 
explains the puzzle of the power -corrections. A namely, the contributions ~ to DIS structure functions, 

usually attributed to higher twists, are known to be present at large but enigmatically disappear at small 
where theycould be extremely impactive. A simple and natural solution of this puzzle was found in Ref. [l6j (see also 
overview []J). In brief, it can be reduced to the following: Among different contributions to the structure functions, 
there are contributions 

T.^^\n-{{Q^+^l'')|^?). (43) 

In the region of large Q^, where > /i^, such contributions can be expanded into power series in the following 
way: 

r„^ln"(QVM^)+5:(^)\ (44) 

The logarithmic term in Eq. (|44l) is included into the leading twist contributions whereas the power terms are 
identical to the power terms attributed to higher twists. Such terms are conventionally supposed to have non- 
perturbative origin. However, Eq. (|44l) explicitly demonstrates that there are the power terms of purely perturbative 
nature. Surely, they should be accounted for in the first time and only the rest should be attributed to higher twists. 

On the other hand, the expansion pij) holds at large only. In the small-Q^ region the power expansion of 
Eq. (03]) looks different: 

r,.~i + i:(2!)\ (45) 

The first term in the r.h.s. of Eq. (|45|) is included into the leading twist contributions while the other terms are again 
the power-Q^ corrections, though with in the nominator. Eq. (|45|) proves that the power corrections ^ 1/(Q^)" 
can never become singular at small Q^. Comparing Eq. (H^ to (H5t reveals that the corrections ^ 1/(Q^)" are about 
vanishing at 

(46) 

where none of expansions (I44I45P can be used. In Ref. [Ts*! we fixed /i for the non-singlet structure functions Ff^ and 
gi^: for them fj,Ns ~ 1 GeV. Our prediction of vanishing power corrections ^ 1/(Q^)" at values of approaching 
~ 1 Gev^ perfectly agrees with results obtained by analysis of experimental data. 

VIII. CONCLUSION 

In the present paper we have analyzed the -dependence of the structure functions, considering separately this 
dependence in the kinematic region of large and small x. Our consideration embraced both large and small values of 
Q^. We argued against the conventional point of view that the Q^-dependence at any and x is originated by the 
DGLAP -parametrization ag — as{Q^) leading to Q^- dependence of the intercepts of the structure functions. We 
showed that this parametrization holds exceptionally in the region of large x and Q^, providing that in this region 
the factorization scale should be large. Outside this region, as should be replaced by the effective coupling cte//- At 
small X, coupling ckg// depends on both transverse and longitudinal momenta, which destroys the factorization of the 
phase space into the longitudinal and transverse spaces taking place in DGLAP. In addition, the integration limits 
do not involve Q^. As a result, the effective coupling does not depend on in the small- a; region. This leads to the 
Regge small- a; asymptotics of the structure functions, with intercepts independent of Q^. We also argued against the 
shift as{Q^) — as{Q^ + fJ.'^) used in the literature to keep as in perturbative domain at small Q^. In this regard, we 
reminded that the parametrization as{Q^) as well as DGLAP in whole should be used at large only. Superficial 
use of the shift as{Q^) as{Q^ + M^) small has led to various misconceptions known in the literature. In 
contrast, we advocated the use of shift + /i^ in such evolution equations, derivation of which is not based 

on assuming large. We stressed that the shift did not involve parametrization of the QCD coupling. Eventually, 
we have considered theoretical grounds for such a shift and discussed its application to the power Q^- contributions 
to the structure functions. 
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